High-accuracy IMEP computational technique using a low-resolution encoder and an indirect integration process

ABSTRACT

A method for computing indicated mean effective pressure (IMEP) in an internal combustion engine using sparse input data. The method uses an indirect integration approach, and requires significantly lower resolution crankshaft position and cylinder pressure input data than existing IMEP computation methods, while providing calculated IMEP output results which are very accurate in comparison to values computed by existing methods. By using sparse input data, the indirect integration method offers cost reduction opportunities for a manufacturer of vehicles, engines, and/or electronic control units, through the use of lower cost sensors and the consumption of less computing resources for data processing and storage.

BACKGROUND OF THE INVENTION

1. Field of the Invention

This invention relates generally to a method for computing meaneffective pressure in an engine and, more particularly, to a method forcomputing indicated mean effective pressure (IMEP) in an internalcombustion engine using an indirect integration method which provides ahighly accurate result even when using a low-resolution crankshaftposition encoder and using less frequent measurement of cylinderpressure input data than required by existing IMEP calculation methods.

2. Discussion of the Related Art

Most modern internal combustion engines employ a number of sophisticatedcontrol strategies to optimize performance, fuel economy, emissions, andother factors. Among the many parameters used to control an engine'soperation, indicated mean effective pressure (IMEP) is one of the moreimportant. IMEP is used as a measure of the amount of work an engine isperforming, or as a measure of the torque that is being provided by theengine. Engine control strategies are often designed around IMEP, and ofcourse these strategies will be effective for controlling the engineonly if IMEP is calculated with a sufficient degree of accuracy.

While methods for calculation of IMEP are known in the art, existingmethods require a high-resolution crankshaft position encoder andfrequent measurement of cylinder pressure data in order to obtain anaccurate IMEP calculation. Requiring high-resolution crankshaft positionand cylinder pressure data has a number of disadvantages, including thecost of the crank position encoder, the cost associated with the digitalmemory required to store the high-resolution cylinder pressure data overtime, and the cost associated with computing power needed in electroniccontrol units in order to process the large amounts of crank positionand cylinder pressure data for IMEP calculations.

A need exists for a method of calculating indicated mean effectivepressure which provides the accuracy needed for proper control of theengine, but which does not require high-resolution crankshaft positionand cylinder pressure data as input. Such a method can provide asignificant benefit in terms of cost savings and simplification for amanufacturer of engines or vehicles.

SUMMARY OF THE INVENTION

In accordance with the teachings of the present invention, an indirectintegration method is disclosed for computing indicated mean effectivepressure (IMEP) in an internal combustion engine using sparse inputdata. The indirect integration method requires significantly lowerresolution crankshaft position and cylinder pressure input data thanexisting IMEP computation methods, while providing calculated IMEPoutput results which are very accurate in comparison to values computedby existing methods. By using sparse input data, the indirectintegration method enables the use of a low-resolution crankshaftposition encoder and requires less computing resources for dataprocessing and storage.

Additional features of the present invention will become apparent fromthe following description and appended claims, taken in conjunction withthe accompanying drawings.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a diagram of a multi-cylinder engine in a vehicle, showing theelements involved in computing indicated mean effective pressure;

FIG. 2 is a flow chart diagram of a first method for calculatingindicated mean effective pressure using sparse input data; and

FIG. 3 is a flow chart diagram of a second method for calculatingindicated mean effective pressure using sparse input data.

DETAILED DESCRIPTION OF THE EMBODIMENTS

The following discussion of the embodiments of the invention directed toan indirect integration method for calculating indicated mean effectivepressure in an engine using sparse input data is merely exemplary innature, and is in no way intended to limit the invention or itsapplications or uses. For example, pumping mean effective pressure andnet mean effective pressure can also be calculated using the methods ofthe present invention.

Engines in most modern automobiles use sophisticated electronic controlunits for the purpose of controlling many parameters of engineoperation—including the amount and timing of fuel injection, sparktiming for spark ignition engines, the amount of exhaust gasrecirculation to be used, and boost pressure for turbo-charged orsuper-charged engines. These parameters and others are preciselycontrolled in an effort to optimize engine performance, fuel economy,and emissions. In many engine controllers, indicated mean effectivepressure (IMEP) is used as an important input parameter to the controlstrategy. IMEP may be thought of as the average pressure over a powercycle in the combustion chamber of an engine, and it is therefore alsorepresentative of the work done by the engine during one cycle, or thetorque being output by the engine over one cycle. IMEP is normallymeasured only during the power cycle of an engine, which includes thecompression stroke and the expansion or power stroke. IMEP can be usedto design an engine control strategy which strives to match the actualtorque being delivered by the engine with the torque being requested bythe driver by way of accelerator pedal position.

Other pressure-related parameters can also be useful in engine controlstrategies. Pumping mean effective pressure (PMEP) is the averagepressure over a pumping cycle (exhaust and intake strokes) in thecombustion chamber of an engine. Net mean effective pressure (NMEP) isthe average pressure over a complete four-stroke cycle in the combustionchamber of an engine. That is, net mean effective pressure is the sum ofindicated mean effective pressure and pumping mean effective pressure.The ensuing discussion and equations are written in terms of IMEP.However, it will be recognized by one skilled in the art that themethods of the present invention are applicable to any calculation ofmean effective pressure (IMEP, PMEP, or NMEP) by simply selecting theintegral range appropriate for the cycle.

FIG. 1 is a diagram showing an internal combustion engine 10 in avehicle 12. The engine 10 includes a plurality of pistons 14 connectedto a crankshaft 16. Each piston 14 travels reciprocally through acylinder 18, while the crankshaft 16 provides output torque to performuseful work, such as driving the vehicle's wheels or charging anelectrical system. In order to calculate IMEP, crankshaft position datafrom a crankshaft position encoder 20 is required, along within-cylinder pressure data from a cylinder pressure sensor 22. The engine10 can include a pressure sensor 22 in every cylinder 18, or as few asone or two cylinder pressure sensors 22 in the entire engine 10. Datafrom the crankshaft position encoder 20 and the cylinder pressure sensor22 are collected by an engine controller 24, which also calculates IMEPand manages engine operation.

A standard definition of IMEP is shown in Equation (1).

$\begin{matrix}{{IMEP} = {\frac{1}{V_{cyl}}{\int{P \cdot {\mathbb{d}V}}}}} & (1)\end{matrix}$Where V_(cyl) is cylinder volume, P is cylinder pressure, dV isincremental cylinder volume, and the integral is taken over an enginepower cycle running from a crank position of −π to +π (or Bottom DeadCenter (BDC) through one revolution back to BDC).

Various methods of calculating IMEP during engine operation are known inthe art. One common IMEP calculation method is the trapezoidalapproximation, where the integral of Equation (1) is discretized insmall increments of volume and summed over an engine power cycle. Thetrapezoidal approximation of IMEP is shown in Equation (2).

$\begin{matrix}{{IMEP} \cong {\frac{1}{V_{cyl}}{\sum\limits_{k = \theta_{0}}^{\theta_{f}}{\frac{P_{k + 1} + P_{k}}{2} \cdot ( {V_{k + 1} - V_{k}} )}}}} & (2)\end{matrix}$Where P_(k) and P_(k+1) are successive cylinder pressure measurements,V_(k) and V_(k+1) are cylinder volume measurements corresponding toP_(k) and P_(k+1), and the summation is taken in increments of k from avalue θ₀ to a value θ_(ƒ).

Although the trapezoidal approximation IMEP calculation of Equation (2)is widely used, it is very sensitive to sampling resolution. That is,the trapezoidal approximation only yields an accurate value of IMEP ifthe pressure and volume increments k are very small—typically 1 degreeof crank rotation or less. The need for high-resolution crankshaftposition and cylinder pressure data means that the crankshaft positionencoder 20 must have high-resolution capability, and it means thatcylinder pressure data must be taken and processed very frequently.While these capabilities exist in engines today, they drive higher costsin the form of the encoder 20 itself, and analog-to-digital conversion,data processing and storage requirements for the volume of cylinderpressure data.

The goal of the present invention is to relax the requirement forhigh-resolution crank position and cylinder pressure data by providing amethod of computing IMEP which is accurate even when the crank positionand cylinder pressure data is measured far less frequently than everydegree of crank angle. This would allow the crankshaft position encoder20 to be a lower-cost, lower-resolution model, and would requiresignificantly less cylinder pressure data to be processed and stored.This in turn would allow the total cost of a pressure-based controlsystem for the engine 10 to be reduced.

In a first embodiment of the present invention, an indirect integrationmethod of computing IMEP in an engine is provided. The indirectintegration method begins with the introduction of a term PV^(n), whereP is pressure, V is volume, and n is the ratio of specific heats. Bydefinition,d(PV ^(n))=V ^(n) dP+nV ^(n−1) PdV  (3)andd(PV)=VdP+PdV  (4)

Rearranging and integrating Equations (3) and (4) yields;

$\begin{matrix}{{\int{P{\mathbb{d}V}}} = {\frac{1}{n - 1} \cdot \lbrack {{\int{\frac{1}{V^{n - 1}} \cdot {\mathbb{d}( {P\; V^{n}} )}}} - {\int{\mathbb{d}( {P\; V} )}}} \rbrack}} & (5)\end{matrix}$

If the integral of Equation (5) is taken over a crank angle range fromθ₀ to θ_(ƒ), Equation (5) can be discretized and written as;

$\begin{matrix}{{\int_{\theta_{0}}^{\theta_{f}}{P{\mathbb{d}V}}} \cong {\frac{1}{n - 1} \cdot \lbrack {{{\frac{1}{{\overset{\_}{V}}^{n - 1}} \cdot {\mathbb{d}( {P\; V^{n}} )}}{\quad }_{\theta_{0}}^{\theta_{f}}} - {( {P\; V} ){\quad }_{\theta_{0}}^{\theta_{f}}}} \rbrack}} & (6)\end{matrix}$

It can be seen that the left-hand side of Equation (6) is the definitionof IMEP from Equation (1), with the exception that the (1/V_(cyl))factor is missing. It therefore follows that IMEP can be approximated asthe right-hand side of Equation (6), multiplied by the (1/V_(cyl))factor, as follows;

$\begin{matrix}{{IMEP} \cong {\frac{1}{V_{cyl} \cdot ( {n - 1} )} \cdot \lbrack {{{\frac{1}{{\overset{\_}{V}}^{n - 1}} \cdot {\mathbb{d}( {P\; V^{n}} )}}{\quad }_{\theta_{0}}^{\theta_{f}}} - {( {P\; V} ){\quad }_{\theta_{0}}^{\theta_{f}}}} \rbrack}} & (7)\end{matrix}$

Equation (7) can then be expanded and written as a summation of discretemeasurements, as follows;

$\begin{matrix}{{IMEP} \cong {\frac{1}{V_{cyl} \cdot ( {n - 1} )} \cdot {\sum\lbrack {{P_{k + \Delta}( {\frac{V_{k + \Delta}^{n}}{{\overset{\_}{V}}_{k}^{n - 1}} - V_{k + \Delta}} )} - {P_{k}( {\frac{V_{k}^{n}}{{\overset{\_}{V}}_{k}^{n - 1}} - V_{k}} )}} \rbrack}}} & (8)\end{matrix}$Where k is the sampling event number, Δ is the increment of crank anglebetween samples, and the remaining terms are as defined above.

It is then possible to define variables G_(k) and H_(k) to represent thevolume terms of Equation (8), as follows;

$\begin{matrix}{G_{k} = {\frac{V_{k + \Delta}}{( {n - 1} )} \cdot ( {\frac{V_{k + \Delta}^{n - 1}}{\frac{1}{\Delta}( {\sum\limits_{j \in {\lbrack{k,{k + \Delta}}\rbrack}}V_{j}^{n - 1}} )} - 1} )}} & (9) \\{H_{k} = {\frac{V_{k}}{( {n - 1} )} \cdot ( {\frac{V_{k}^{n - 1}}{\frac{1}{\Delta}( {\sum\limits_{j \in {\lbrack{k,{k + \Delta}}\rbrack}}V_{j}^{n - 1}} )} - 1} )}} & (10)\end{matrix}$

It is notable that the variables G_(k) and H_(k) contain only constantsand volume-related terms, which are known functions of cylinder volumeand crank position. Therefore G_(k) and H_(k) can be computed offlineand stored for any particular engine geometry, as they do not depend oncylinder pressure or any other real-time engine performance factor.

Substituting G_(k) and H_(k) into Equation (8) yields;

$\begin{matrix}{{IMEP} \cong {\frac{1}{V_{cyl}}{\sum\lbrack {{P_{k + \Delta} \cdot G_{k}} - {P_{k} \cdot H_{k}}} \rbrack}}} & (11)\end{matrix}$

Again, it is noteworthy that V_(cyl) is a constant, and the terms G_(k)and H_(k) are pre-computed and known for each sampling event k.Therefore, IMEP can be calculated using Equation (11) by simplymultiplying a cylinder pressure measurement, P_(k+Δ), by itsvolume-related term, G_(k), subtracting the product of the previouscylinder pressure measurement, P_(k), and its volume-related term,H_(k), and summing the results over an engine power cycle.

FIG. 2 is a flow chart diagram 40 of the indirect integration method forcomputing IMEP discussed in the preceding paragraphs. At box 42, initialvalues are defined, where Δ is the sampling resolution, n is set equalto 1.4 which is the normal specific heat ratio for air, and values ofV_(k) are calculated as the cylinder volume as a function of the crankangle θ at each sampling event k. At box 44, values for thevolume-related terms G_(k) and H_(k) are computed as a function of thecrank angle θ at each sampling event k. The calculations of the box 44are also completed in an initialization phase prior to real-time engineoperation, as the calculations are a function only of engine geometryand the chosen crank angle increment Δ. During engine operation, thereal-time calculation of IMEP is handled at box 46 using Equation (11)in a summation over one engine power cycle, where the cylinder pressuredata is sampled at each crank angle θ corresponding to a crank angleincrement Δ, as shown at box 48. The pressure data in the box 48 ismeasured by cylinder pressure sensors 22. At the end of one completeengine power cycle, the value of IMEP is output at box 50 as the resultof the summation at the box 46. The IMEP value from the box 50 is thenused by the engine controller 24 to control engine operation, asdiscussed previously.

In a second embodiment of the present invention, a cubic splineintegration method of computing IMEP in an engine is provided. In thecubic spline integration method, a cubic spline is fitted to theintegral Equation (1). This allows IMEP to be calculated with sufficientaccuracy, even when using sparse cylinder pressure data. According tothis method, ƒ(x) is defined as a continuous function, as follows;

$\begin{matrix}{{f(x)} = {\frac{1}{V_{cyl}} \cdot P \cdot \frac{\mathbb{d}V}{\mathbb{d}\theta}}} & (12)\end{matrix}$Where V_(cyl) is cylinder volume, P is cylinder pressure, and dV/dθ isthe first derivative of cylinder volume with respect to crank angleposition θ.

The function ƒ(x) is defined to have a continuous third derivativethrough the interval [a, b], where;a=x ₀ <x ₁ < . . . <x _(n−1) <x _(n) =b  (13)

It can be seen from Equation (1) and Equation (12) that a value for IMEPcan be obtained by integrating the function ƒ(x) over one power cycle,that is, from;

$\begin{matrix}{{x_{0} = {{- 180}{{^\circ} \cdot \frac{\pi}{180{^\circ}}}}}{{to};}} & (14) \\{x_{n} \cong {{+ 179}{{^\circ} \cdot \frac{\pi}{180{^\circ}}}}} & (15)\end{matrix}$

Therefore an equation for IMEP can be written as;

$\begin{matrix}{{IMEP} = {S_{\theta_{f}} \cong {\frac{1}{V_{cyl}}{\int_{\theta_{0}}^{\theta_{f}}{P\ {\mathbb{d}V}}}}}} & (16)\end{matrix}$Where the function S is the cubic spline integral of ƒ, θ₀=x₀ andθ_(ƒ)=x_(n).

An algorithm for computing IMEP via the cubic spline integral S isdefined as follows. First, a function M is defined as the firstderivative of ƒ. Solving for M at the initial point θ₀ yields;

$\begin{matrix}{M_{0} = {{f^{\prime}( \theta_{0} )} = {\frac{1}{V_{cyl}} \cdot ( {{{P^{\prime}( \theta_{0} )}\frac{\mathbb{d}V}{\mathbb{d}\theta}( \theta_{0} )} + {{P( \theta_{0} )}\frac{\mathbb{d}^{2}V}{\mathbb{d}^{2}\theta}( \theta_{0} )}} )}}} & (17)\end{matrix}$

In Equation (17), θ₀ is the beginning of the power cycle, that is, thebeginning of the compression stroke, which is at a crank position ofBottom Dead Center (BDC), or −π. At this point, the cylinder pressurecan be approximated as constant, and therefore;P′(θ₀)=0  (18)P(θ₀) can be easily obtained from the cylinder pressure sensor 22.

In order to resolve the first and second derivative terms of Equation(17), a formulation for volume V as a function of crank angle θ isneeded. This can be expressed as follows;V(θ)=K₁−K₂(cos(θ)+√{square root over (R²−sin²(θ)))}  (19)Where K₁ and K₂ are engine-related constants, and R is defined as r/L,with r being the crank radius and L being the connecting rod length.From Equation (19), the calculation of dV/dθ and d²V/d²θ becomestraightforward to one skilled in the art.

Next, a crank angle increment h is defined such that;h=θ _(i)−θ_(i−1)  (20)Where h can be defined as any value which may be suitable for thepurpose, and i is step number. Since the objective of this method is tocompute a value of IMEP using sparse cylinder pressure data, values of hwhich are significantly larger than 1 degree of crank angle will beexplored, such as 3 degrees or 6 degrees.

Now a recursive calculation can be set up, where each power cycle of theengine 10 begins by initializing;

$\begin{matrix}{M_{0} = {{f^{\prime}( \theta_{0} )} = {\frac{1}{V_{cyl}} \cdot ( {{P( \theta_{0} )}\frac{\mathbb{d}^{2}V}{\mathbb{d}^{2}\theta}( \theta_{0} )} )}}} & (21)\end{matrix}$Where P(θ₀) is the measured cylinder pressure at the cycle initiationlocation of Bottom Dead Center, V_(cyl) is total cylinder volume, and(d²V/d²θ)(θ₀) is the second derivative of Equation (19) evaluated at thecycle initiation location of BDC. Also, at the cycle initiation, ƒ₀=0because the factor dV/dθ is zero at BDC, and S₀=0 by definition.

Then, for each step i of crank angle increment h, the functions ƒ and Mcan be solved sequentially as follows;

$\begin{matrix}{{f_{i} = {{f( \theta_{i} )} = {{\frac{1}{V_{cyl}} \cdot {P( \theta_{i} )} \cdot \frac{\mathbb{d}V}{\mathbb{d}\theta}}( \theta_{i} )}}}{{and};}} & (22) \\{M_{i} = {\frac{2( {f_{i} - f_{i - 1}} )}{h} - M_{i - 1}}} & (23)\end{matrix}$Where P(θ_(i)) is the measured cylinder pressure at the current step i,(dV/dθ)(θ_(i)) is the first derivative of V with respect to θ evaluatedat the current step i, and h is the crank angle increment.

Then the cumulative cubic spline function S can be calculated from theprevious value of S, the current and previous values of M, and theprevious value of ƒ, as follows;

$\begin{matrix}{S_{i} = {S_{i - 1} + {\frac{1}{6}{h^{2}( {M_{i} + {2M_{i - 1}}} )}} + {hf}_{i - 1}}} & (24)\end{matrix}$The function S is calculated in a cumulative fashion from a value ofS₀=0 at cycle initiation until the power cycle ends when θ_(i)=θ_(ƒ),which is at BDC at the end of the power stroke. At that point, IMEP forthe completed power cycle is output as the final value of S; that is,IMEP=S_(θ) _(ƒ) . Then a new cycle is initiated.

FIG. 3 is a flow chart diagram 80 of the cubic spline integration methodof computing IMEP discussed in the preceding paragraphs. One-timeinitialization calculations are handled at box 82. The values computedat the box 82 are constants associated with a particular engine design,such as stroke, connecting rod length, piston area, and cylinder volume.At box 84, cycle initiation calculations take place. These calculationsinclude measuring cylinder pressure P and calculating the first andsecond derivatives of V and the functions ƒ, M, and S—all at the cycleinitiation location of θ=θ₀=−π, which is at BDC before the compressionstroke. At box 86, cylinder pressure P is measured, the first derivativeof V is calculated, and the function ƒ is evaluated for each step i. Atbox 88, the functions M and S are evaluated per Equations (23) and (24).At decision diamond 90, the crank angle θ is checked to see if the powercycle has been completed. If θ_(i)≧θ_(ƒ), then at box 92 the value ofIMEP for the cycle is output as the final value of S, and a new cycle isstarted at the box 84. If θ_(i)<θ_(ƒ) at the decision diamond 90, thenthe current cycle calculations continue at the box 86 for the next stepi.

Both the indirect integration method and the cubic spline integrationmethod of computing IMEP have been tested with simulations using realengine data. IMEP calculations using the disclosed methods with sparsedata (sampling resolution at crank rotation increments of 3 degrees and6 degrees) were found to be within 2% of IMEP calculations using densedata (crank rotation increment of 1 degree) in a traditional trapezoidalapproximation. This variance of less than 2% is well within anacceptable range for using IMEP in the engine controller 24. Evensampling resolutions as large as 10 degrees were found to produceacceptable IMEP results using the disclosed methods. By using cylinderpressure data at crank position increments of 6 degrees instead of 1degree, the disclosed methods achieve the desired goal of relaxing therequirements of high-resolution crank position and cylinder pressuredata, and enable a reduction in the total cost of pressure-based enginecontrol systems.

The foregoing discussion discloses and describes merely exemplaryembodiments of the present invention. One skilled in the art willreadily recognize from such discussion and from the accompanyingdrawings and claims that various changes, modifications and variationscan be made therein without departing from the spirit and scope of theinvention as defined in the following claims.

What is claimed is:
 1. A method for computing mean effective pressure inan internal combustion engine, said method comprising: defining asampling resolution for a series of sampling events, said samplingresolution being an amount of crankshaft rotation between the samplingevents; computing and storing a volume array, said volume arraycontaining combustion chamber volume as a function of crankshaftposition for each sampling event; computing and storing arrays of afirst and a second volume-related variable, said volume-relatedvariables including values for each crankshaft position corresponding tothe sampling resolution; where the first volume-related variable isdefined by the equation:$G_{k} = {\frac{V_{k + \Delta}}{( {n - 1} )} \cdot ( {\frac{V_{k + \Delta}^{n - 1}}{\frac{1}{\Delta}( {\sum\limits_{j \in {\lbrack{k,{k + \Delta}}\rbrack}}V_{j}^{n - 1}} )} - 1} )}$and the second volume-related variable is defined by the equation:$H_{k} = {\frac{V_{k}}{( {n - 1} )} \cdot ( {\frac{V_{k}^{n - 1}}{\frac{1}{\Delta}( {\sum\limits_{j \in {\lbrack{k,{k + \Delta}}\rbrack}}V_{j}^{n - 1}} )} - 1} )}$where k is sampling event number, n is a ratio of specific heats, Δ isthe sampling resolution, and V is combustion chamber volume from thevolume array evaluated at any sampling event k; running the engine andtaking a cylinder pressure measurement at each crankshaft positioncorresponding to the sampling resolution; storing the cylinder pressuremeasurement for a current sampling event and the cylinder pressuremeasurement for a previous sampling event for calculation purposes; andcomputing mean effective pressure for an engine cycle using an indirectintegration calculation, where the indirect integration calculation is afunction of the cylinder pressure measurements, the first volume-relatedvariable, and the second volume-related variable.
 2. The method of claim1 wherein the indirect integration calculation of mean effectivepressure is defined by the equation:${IMEP} = {\frac{1}{V_{cyl}}{\sum\lbrack {{P_{k + \Delta} \cdot G_{k}} - {P_{k} \cdot H_{k}}} \rbrack}}$where k is the sampling event number, Δ is the sampling resolution,V_(cyl) is cylinder volume of a cylinder in the engine, P_(k+Δ) is thecylinder pressure measurement at the current sampling event, G_(k) isthe first volume-related variable evaluated at the previous samplingevent, P_(k) is the cylinder pressure measurement at the previoussampling event, H_(k) is the second volume-related variable evaluated atthe previous sampling event, and the summation is calculated over acycle of the engine.
 3. The method of claim 1 wherein the ratio ofspecific heats is 1.4.
 4. The method of claim 1 wherein the samplingresolution is three degrees or greater of crankshaft rotation.
 5. Themethod of claim 1 wherein the sampling resolution is six degrees orgreater of crankshaft rotation.
 6. The method of claim 1 furthercomprising using the calculated value of mean effective pressure in anengine controller to control operation of the engine.
 7. A method forcomputing and using indicated mean effective pressure in an internalcombustion engine, said method comprising: defining a samplingresolution for a series of sampling events, said sampling resolutionbeing an amount of crankshaft rotation between the sampling events;computing and storing a volume array, said volume array containingcombustion chamber volume as a function of crankshaft position for eachsampling event; computing and storing arrays of a first and a secondvolume-related variable, said volume-related variables including valuesfor each crankshaft position corresponding to the sampling resolution;where the first volume-related variable is defined by the equation:$G_{k} = {\frac{V_{k + \Delta}}{( {n - 1} )} \cdot ( {\frac{V_{k + \Delta}^{n - 1}}{\frac{1}{\Delta}( {\sum\limits_{j \in {\lbrack{k,{k + \Delta}}\rbrack}}V_{j}^{n - 1}} )} - 1} )}$and the second volume-related variable is defined by the equation:$H_{k} = {\frac{V_{k}}{( {n - 1} )} \cdot ( {\frac{V_{k}^{n - 1}}{\frac{1}{\Delta}( {\sum\limits_{j \in {\lbrack{k,{k + \Delta}}\rbrack}}V_{j}^{n - 1}} )} - 1} )}$where k is sampling event number, n is a ratio of specific heats, Δ isthe sampling resolution, and V is combustion chamber volume from thevolume array evaluated at any sampling event k; running the engine andtaking a cylinder pressure measurement at each crankshaft positioncorresponding to the sampling resolution; storing the cylinder pressuremeasurement for a current sampling event and the cylinder pressuremeasurement for a previous sampling event for calculation purposes;computing indicated mean effective pressure for an engine cycle using anindirect integration calculation, where the indirect integrationcalculation is a function of the cylinder pressure measurements, thefirst volume-related variable, and the second volume-related variable;and using the calculated value of indicated mean effective pressure inan engine controller to control operation of the engine, includingcontrol of fuel flow to the engine.
 8. The method of claim 7 wherein theindirect integration calculation of indicated mean effective pressure isdefined by the equation:${IMEP} = {\frac{1}{V_{cyl}}{\sum\lbrack {{P_{k + \Delta} \cdot G_{k}} - {P_{k} \cdot H_{k}}} \rbrack}}$where k is the sampling event number, Δ is the sampling resolution,V_(cyl) is cylinder volume of the engine, P_(k+Δ) is the cylinderpressure measurement at the current sampling event, G_(k) is the firstvolume-related variable evaluated at the previous sampling event, P_(k)is the cylinder pressure measurement at the previous sampling event,H_(k) is the second volume-related variable evaluated at the previoussampling event, and the summation is calculated over a complete powercycle of the engine.
 9. The method of claim 8 wherein the ratio ofspecific heats is 1.4.
 10. The method of claim 9 wherein the samplingresolution is three degrees or greater of crankshaft rotation.
 11. Themethod of claim 10 wherein the sampling resolution is six degrees orgreater of crankshaft rotation.
 12. A system for computing and usingindicated mean effective pressure in an internal combustion engine, saidsystem comprising: a low-resolution crankshaft position encoder formeasuring crankshaft position at a sampling resolution, where thesampling resolution is greater than one degree of crankshaft rotation; acylinder pressure sensor for measuring cylinder pressure at eachcrankshaft position corresponding to the sampling resolution; and anengine controller configured to collect data from the crankshaftposition encoder and the cylinder pressure sensor, calculate indicatedmean effective pressure using an indirect integration algorithm, and usethe calculated indicated mean effective pressure to control operation ofthe engine, where the indirect integration algorithm includes defining asampling resolution, said sampling resolution being the amount ofcrankshaft rotation between sampling events, computing and storing avolume array, said volume array containing combustion chamber volume asa function of crankshaft position for each sampling event, computing andstoring arrays of a first and a second volume-related variable, saidvolume-related variables including values for each crankshaft positioncorresponding to the sampling resolution, running the engine and takinga cylinder pressure measurement at each crankshaft positioncorresponding to the sampling resolution, storing the cylinder pressuremeasurement for a current sampling event and the cylinder pressuremeasurement for a previous sampling event for calculation purposes, andcalculating indicated mean effective pressure for an engine cycle as afunction of the cylinder pressure measurements, the first volume-relatedvariable, and the second volume-related variable.
 13. The system ofclaim 12 wherein the first volume-related variable is defined by theequation:$G_{k} = {\frac{V_{k + \Delta}}{( {n - 1} )} \cdot ( {\frac{V_{k + \Delta}^{n - 1}}{\frac{1}{\Delta}( {\sum\limits_{j \in {\lbrack{k,{k + \Delta}}\rbrack}}V_{j}^{n - 1}} )} - 1} )}$and the second volume-related variable is defined by the equation:$H_{k} = {\frac{V_{k}}{( {n - 1} )} \cdot ( {\frac{V_{k}^{n - 1}}{\frac{1}{\Delta}( {\sum\limits_{j \in {\lbrack{k,{k + \Delta}}\rbrack}}V_{j}^{n - 1}} )} - 1} )}$where k is sampling event number, n is a ratio of specific heats, Δ isthe sampling resolution, and V is combustion chamber volume from thevolume array evaluated at any sampling event k.
 14. The system of claim13 wherein the indicated mean effective pressure is calculated using theequation:${IMEP} = {\frac{1}{V_{cyl}}{\sum\lbrack {{P_{k + \Delta} \cdot G_{k}} - {P_{k} \cdot H_{k}}} \rbrack}}$where k is the sampling event number, Δ is the sampling resolution,V_(cyl) is cylinder volume of the engine, P_(k+Δ) is the cylinderpressure measurement at the current sampling event, G_(k) is the firstvolume-related variable evaluated at the previous sampling event, P_(k)is the cylinder pressure measurement at the previous sampling event,H_(k) is the second volume-related variable evaluated at the previoussampling event, and the summation is calculated over a complete powercycle of the engine.
 15. The system of claim 13 wherein the ratio ofspecific heats is 1.4.
 16. The system of claim 12 wherein the samplingresolution is three degrees or greater of crankshaft rotation.
 17. Thesystem of claim 12 wherein the engine controller uses the calculatedindicated mean effective pressure to control fuel flow to the engine.